function [p q] = TimeDiscretizers(dt,option)
% Following Algorithms are implemented in the current study
% 1. Fully Implicit Theta Scheme
% 2. Fully Implicit BDF Two parameter Scheme
% 3. Explicit Leap Frog Scheme
% 4. Semi Implicit BDF - By adjusting parameters = Semi Implicit
% Trapezoidal
% 5. Semi Implicit Trapezoidal with Upwinding Discretization
% 6. Semi Implicit Trapezoidal with Central Discretization
% 7. Fully Implicit with Upwinding
global Ux Uy  h  g  p  q r tmax Lx Ly dx dy L Nx Ny y nt;

%Reinitialize the variables

pi1=1;
qi=1;
%r=zeros(Nx*Ny,1);
for j=1:Nx
    for i=1:Ny
        index = (j-1)*Ny+i;
        p(index,1) = pi1*cos(4*pi*(j-1)*dx/Lx);
        q(index,1) = qi;
       % r(index,1)=0;
    end
end


nt = round(tmax/dt);   % number of time steps
%h1 = plot(xplot,p,'yDataSource','p');


if (option ==1)
    disp('1. Chorin Projection Scheme\n');
    Chorin_Projection(dt);

end




